Quantum algorithm for measuring the energy of n qubits 
with unknown pair-interactions 
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The well-known algorithm for quantum phase estimation requires that the considered unitary is 
available as a conditional transformation depending on the quantum state of an ancilla register. We 
present an algorithm converting an unknown n-qubit pair-interaction Hamiltonian into a conditional 
one such that standard phase estimation can be applied to measure the energy. Our essential 
assumption is that the considered system can be brought into interaction with a quantum computer. 
For large n the algorithm could still be applicable for estimating the density of energy states and 
might therefore be useful for finding energy gaps in solid states. 



I. INTRODUCTION 

Finding the energy spectrum of a given Hamiltonian 
is an important task in physics since it determines the 
dynamical and thermodynamical behavior of a quantum 
system. In solid-states physics, for instance, the spec- 
trum, in particular its gaps, are relevant for the transport 
properties and in the famous BCS-theory ]jj the energy 
gap between the ground state and the first excited state 
is decisive for understanding super-conductivity. 

Therefore a lot of efforts has been done on calculat- 
ing the spectrum of many-particle Hamiltonians. Some 
models are known, where the Hamiltonian can explicitly 
be diagonalized Q and advanced mathematical tools like 
non-commutative geometry could be shown to be useful 
for finding energy gaps in perfect crystals or those with 
defects and quasi-periodicity ||. Nevertheless, the diago- 
nalization of generic many-particle Hamiltonians is com- 
putationally hard and even restricted information about 
its spectrum is difficult to obtain, since the dimension of 
the Hilbert space grows exponentially with the number 
of particles. Here we discuss how to obtain information 
about the spectrum of a Hamiltonian by a quantum com- 
puter. Apart from the fact that the algorithm is able to 
find energy gaps in many-particle systems efficiently, it 
has a property which is impossible in principle without 
processing quantum information: we can gain informa- 
tion about the spectrum of a physical system with an 
unknown Hamiltonian provided that there is an inter- 
face between the quantum register and the considered 
system for exchanging quantum information^]. If such 
an interface does not exist the algorithm can be used if 
the Hamiltonian is known and the quantum computer is 
able to simulate the corresponding unitary evolution, i.e., 
if efficient implementations of exp(—iHt) for each t > 
exist. Then the algorithm can use those implementations 



as black-box subroutines^. 

We restrict our attention to an unknown Hamiltonian 
H of the following form. Let H act on the Hilbert space 
H := (C 2 )®", which we shall refer to as the target regis- 
ter. We assume H to consist of 1-qubit terms and pair- 
interactions between the n two-level systems, i.e. 



H := 



j<n,a 



E 

k<l<n,ct,/3 



j k I 



(1) 



where a = x,y, z and o\ is the Pauli matrix a a acting 
on qubit k. The in x in- matrix J and the 3n- vector r 
specify the interaction uniquely 

First we should discuss why the problem can not be 
solved by a simple application of the well-known quan- 
tum algorithm for phase estimation |j],[lO| . The reason is 
that the algorithm for estimating eigenvalues of a unitary 
u does not work with black-box queries of u, it relies on 
the quantum transformation 

u \=u® |1)(1| + 1® |0)(0| 

implementing u if and only if an ancilla qubit is in the 
state |1), i.e., the ancilla qubit controls the implementa- 
tion of u. This causes severe problems if u := exp(— iHt) 
is the natural time evolution, since phase estimation 
would require a system with Hamiltonian 

H® |1)(1|, 

i.e., the time evolution should be switched on and off 
by the ancilla's state. First of all it is not clear how to 
'switch off' the natural dynamics of a system, for example 
the interaction between nuclear spins in a real molecule. 
However, this problem can be solved if certain approx- 
imations are allowed as it is the case in standard de- 
coupling techniques in Nuclear Magnetic Resonance . 
Assume that local transformations of the form 
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1 Of course information about the spectrum of a system's Hamiltonian can be obtained by sophisticated experiments in principle 
without using our algorithm. But performing sophisticated quantum measurements is quantum information processing. 
2 The problem of simulating pair-interaction Hamiltonians by a given one has recently been discussed in Jtl pf . 
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v := v\ ® V2 ® • • • ® v n (2) 

can be applied arbitrarily fast ('fast control limit'). If 
v is applied before the system evolves according to its 
Hamiltonian H and is applied afterwards, then the 
system evolves as if it was subjected to the conjugated 
Hamiltonian vHvr . Concatenations of such 'conjugated 
evolutions' for small time intervals make the system ap- 
proximatively evolving as if it was subjected to the 'aver- 
age Hamiltonian' which is given by the convex combina- 
tion of the different conjugated Hamiltonians Hf7|. This 
technique can be applied for 'switching off' Hamiltonians 
. At first sight it might seem as if this technique could 
be applied to control H by the ancilla's state. But this 
would require to substitute the sequence of decoupling 
local transformations by conditional transformations of 
the form 

v := v® |0)(0| + 1 <g> |1)(1|. 

with v as eq. @). The assumption that v could be im- 
plemented arbitrarily fast is much less justified than the 
usual fast control limit, since such a transformation refers 
to an interaction between ancilla and the target qubits 
and we have no reason to assume that this interaction is 
considerably stronger than the interaction which we want 
to switch off. 

In other words, there is no obvious way to convert H 
into H 1 1) Below we present a scheme converting H 
into H a z , which has essentially the same effect. Our 
assumptions are the following. 

1. The target register Tt can be brought into interac- 
tion with a quantum computers register TZ. This 
interaction acts on the same time scale as H . The 
interaction between Ti and TZ can be switched off 
without disturbing H . 

2. During that time period where TL interacts with 
TZ the Hamiltonian evolution H can be switched 
off. This decoupling is controlled by classical sig- 
nals and not by a state of a quantum register. 

In order to illustrate this assumptions we consider 
a molecule with k := n + m nuclear spins. Assume 
that the interactions between the spin pairs (J, I) for 
j, I < n are unknown and the interactions between all 
the other pairs are known. Given an arbitrary subset 
AI C {1,2,. . . , n + m}, the following procedure can be 
applied for switching off all those terms in the Hamiltoni- 
ans concerning qubits in the set M. For each qubit in M 
choose vectors x\,. . . ,xi of equal dimension d with the 
following property: The entries are 1, x, y, z and for each 
pair Xi, Xj of vectors each pair (a, (3) with a,/3 = l, x, y, z 
appear equally often in the list {x\, Xj) r <d if x\ is the r th 
entry of the vector X4. The vectors x\, . . . ,Xi are said to 
form an orthogonal array If %j = & than the time 

evolution of spin j is conjugated by the unitary trans- 
formation a a during the r th time period, where we have 



used the convention that o\ is the identity map. If the 
time periods are small, the resulting time evolution is 
approximatively the identity for all the spins in M and 
those terms of H which do not involve spins in M are 
unchanged. 

Then we consider the first n spins as the target reg- 
ister with the unknown Hamiltonian H and the other 
spins form the quantum computer's register TZ. The total 
Hamiltonian will be of the same form as eq. ([!]) charac- 
terized by a 3(n + m)-vector r and a 3(n + m) x 3(n + rn)- 
matrix J. 

In order to switch off the interaction between Ti. and TZ 
we choose M := {n + 1, . . . , n + m}. Then the unknown 
Hamiltonian H on the spins 1, . . . , n is remaining. 

If the Hamiltonian on the target register should be 
switched off in such a way that the interaction between 
one specific qubit j < n in the target register and one 
specific qubit n + I in the ancilla register remains, one 
has to take M as the complement of the set {j, n + I}. 
Note that this does not remove the 1-qubit terms <7™ +J 
and o~ 3 a . This is important since o~ J a is unknown and is 
therefore disturbing if we want to implement a definite 
2-qubit transformation on the qubit pair (J, n+l). There- 
fore we have to switch it off. This can be done as follow. 
The interaction between j and n + / is given as 

Ei n+l 
g a ,l30- J a o- p . 

a, /3 

Choose a specific pair a, (3 such that g a .p ^ 0. Cancel 
all the terms cr~ er^ +i with (a, (3) ^ (a, 13) by conjugation 

of the evolution with the 4 unitaries 1, a J a , o~ J a a^ +l 
on 4 time periods of equal length. The result is that the 
term 

rJ rri 4- r n+l rr n+l 4- n ol r^ rr n+l 

is remaining. The 1-qubit terms can be cancelled by con- 
jugating the evolution with cr^,CT^, +i for half of the time 
period with a 1 7^ a and (3 1 ^ (3. Since this changes 
the sign of the operators o~ J a (by the anti-commutation 
property of the Pauli- matrices) and cr^, the bilinear term 
o~ 3 a o-p is unchanged. 

This shows that the assumptions 1. and 2. above are 
justified. Note that it is not relevant which interaction 
between ancilla qubits and register is available. If the 
true physical interaction is 

E i l+n 

a,f3 

then it can be converted into each other Hamiltonian 
with coefficients g a ^ (see Q) in the sense of the 'average 
Hamiltonian' method. Hence, if we write 'switch on the 
interaction o~ x ® o~z we do not assume the real physical 
Hamiltonian to be of this form, we assume only that the 
true Hamiltonian of the system can be used for simulat- 
ing the required term. This is the continuous analogue of 
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the usual way of describing algorithms by basic gates as 
'controlled not' operations: it does not matter whether 
the quantum computer really has the 'controlled not' as a 
basic operation but it should be capable of 'simulating' it 
by those transformations which are really available. Our 
way of describing the algorithm by 'switching on and 
off interactions' should therefore only be considered as a 
convenient language for continuous quantum algorithms. 

II. THE ALGORITHM 

First we consider only one ancilla qubit and describe 
how to convert the Hamiltonian H®1 into the conditional 
Hamiltonian H ® a z controlled by the ancilla's state. 

We can switch off all the coupling except of those on 
one specific qubit pair (j, k) by so-called selective decou- 
pling ||. Furthermore we can select those terms in the 
remaining Hamiltonian which contain only a J a and o~% for 
specific a and (3. The remaining term is given by 

Hj, k ,a,i3 ■= ri<ri + r h a k g + Jj,k, a ,0(T j a (J k B 

First we convert Hj^, a ,/3 into the term 

This can be done by dividing the considered small time 
interval [0, S] into two intervals [0,(5(1/2 — l/(n — 1))] 
and [(5(1/2 — l/(n — 1 )),£]. During the first interval the 
evolution according to Hj t k,a,p is conjugated by & 3 a <Jp. 
This reverses the sign of the 1-qubit terms whereas the 
2-qubit term is unchanged. During the second interval 
the evolution according to Hj^ >a> p is applied. 

Now we can describe the conversion into a controlled 
Hamiltonian. By identifying x, y, z with the additive 
group F-3 := {0, 1, 2} it can be explained as follows. 

1. Apply the unitary transformation uj defined by the 
equation 

Uj(T J a u] = a 3 a+1 . 

2. Apply time evolution corresponding to H j k a g for 
a small time period e. 

3. Switch on the interaction between target and an- 
cilla qubit given by 

H J ■= ® °» 

for the time e. 

4. Simulate evolution according to — ifj k a g for the 
time e. Inverting unknown Hamiltonians is de- 
scribed in ||. In our case the inversion subroutine 
can be reduced as follows. Choose a large number 
p. Divide the time interval [0, 3e] into 3p intervals 



of equal length. During the interval I conjugate 
the evolution according to H / k a g by the unitary 

transformations o- J a+1 , cr^ +1 , o"„ +1 cr^ +1 if I is 0, 1, or 
2 mod 2, respectively. In the limit p — > oo, we ob- 
tain a perfect simulation of the evolution according 

5. Switch on the interaction 

for the time e. 

6. Apply u|. 

7. Repeat steps 1 to 6 with qubit k. 

Up to an error of order e 3 this simulates the time evo- 
lution corresponding to 

iu A H j,Ka,^ Hj}u] + iu k [H' :jXafB ,H j \u ] k = 2Hl ka B <g> a z 
for the time e 2 with the definition 

H",k,a,0 ■= ^—[(« + r k a k p ) + Jj, k , a ^i4- 
Due to the equation 

j</c<n,a,/5 

the concatenation of the above procedure for all un- 
ordered pairs (J, k) and all a, simulates time evolution 
according to 

H <g> a z . 

Note that the time required for obtaining the unitary 
transformation 

exp(—iH ® cr z S) 

goes to zero for (5 — > but the time overhead required for 
this simulation depends on the desired accuracy since the 
resulting evolution is only of second order a non-trivial 
one. Note that cr^crg can also be converted into o- J a ap®a z 
by a conjugation of the evolution with an appropriate 2- 
qubit gate acting on qubit k and the ancilla qubit, since 
o~p (8 1 and a k a ® a z have the same spectrum as operators 
on 2 qubits. But this kind of converting the Hamiltonian 
into a conditional one does not have the property that 
an infinitesimal time step of the algorithm requires only 
infinitesimal time. Whether this is a true disadvantage 
for the desired accuracies has to be checked by thorough 
numerics. 

Now we can apply usual phase estimation procedure 
to the conditional Hamiltonian evolution. 

Now we initialize the ancilla register by a Hadamard 
transformation on each qubit. We implement the trans- 
formation 
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exp(-i(H ® <r J z 2 j T)) 

for each ancilla qubit j — 0, . . . , m — I. The total run- 
ning time of this implementation is (2 m — l)r times the 
time overhead for simulating H ® a z with appropriate 
accuracy. 

We choose r in such a way that 2r A = ir if A is an up- 
per bound for the difference between greatest and small- 
est eigenvalue of H given by prior knowledge. For every 
known systems in many-particle physics, A is always of 
the order of n, since energy per particle is a well-defined 
quantity in the thermodynamic limit n — > oo. The size 
to of the ancilla register has to be chosen in such a way 
that 2~ m r is smaller than the desired accuracy of the 
energy measurement. It is convenient to consider the an- 
cilla qubit j as the j th number of a binary digit. Then 
the initial state of TZ can be written as 

./1m — ' 
VZ Z<2>™ 

If the target register is in an eigenstate of H with eigen- 
value E we obtain (up to an irrelevant global phase) the 
ancilla state 

a/2™ ~ 

since exp(— i<J z r) produces a relative phase difference 
exp(— 2£t) between the ancilla states |0) and |1). 

By inverse Fourier transformation, one obtains the 
state 

±_J2e- i2 " kl ^e i2ET \k). 
k,l 

Measuring the ancilla register in the standard basis pro- 
vides good estimations for E since the probability that 
the measured result k differs more than e from 2 m ET/ir 
is less than l/(2e — 2) ||. If 2 _m /r is much smaller than 
the energy gaps of H, the algorithm projects onto the 
eigenstates of H [ fTo) , For large n, we cannot expect this 
to be achievable since the gaps decrease exponentially 
with n in the generic case. However, in this case the al- 
gorithm can be used for estimating the density of energy 
states, since the measured results will mostly be around 
those energy values where the density is high, provided 
that the target register has been initialized in the max- 
imally mixed state. Estimating the density of states is 
an important task in solid-state physics [l2] ] and spectral 
gaps can be detected by this method. Most interesting 
energy gaps in solid states physics do not depend on the 
number of particles |jl3|,[l) . They could be detected with- 
out increasing the size of the ancilla register. But one 
should emphasize that only those gaps can be detected 



which are localized around typical energy values accord- 
ing to the initial density matrix. Therefore it might be 
more useful to start with thermal equilibrium states in 
the target register in order to find gaps around values 
which are typical for the desired temperature^]. 

We should mention a difficulty which appears for large 
n and how it could be overcome. We assumed that each 
qubit in the target register can be brought into interac- 
tion with each qubit of the ancilla register. Keeping in 
mind the example with nuclear spins in a molecule, one 
should take into account that the strength of the interac- 
tions is strongly decreasing with the distance between the 
spins. Accordingly, the interaction time between ancilla 
and target spin has to be increased for large distances 
and the resulting running time might become unaccept- 
able. Therefore we emphasize that the algorithm can be 
rewritten in such a way, that it is not necessary to have 
an interaction between each qubit in Ti and each qubit 
in TZ. Instead, one can transfer the information qubit by 
qubit into a quantum register. To speak more precisely, 
one can realize the algorithm within the following setting 
instead of assumptions 1) and 2) above: 

1. Apart from TZ, there is a quantum register Ti' with 
n qubits such that the 'information exchange' 

w:\<j>)<& \ip) <-> \ip) <8> \4>) 

can be implemented at definite times in such a way 
that exp(-iHt) can be implemented on Ti' by the 
following subroutine: 

• implement w. 

• wait the time t. 

• implement w. 

2. The system Ti' ®TZ has full capabilities of quantum 
computation. 

Then Ti can be substituted by Ti' if the instruction 
'wait the time V is substituted by the subroutine above. 
Due to the quantum computation capabilities on Ti' ® TZ 
we can clearly implement any desired transformation of 
the form exp(— ia a <g> a z t) for all t € K between any ar- 
bitrary qubit pair. This shows once more that instruc- 
tions of the form 'switch on the interaction (g> cr^ +i ' as 
above should not be taken to literally. Having full capa- 
bilities of quantum computation, we can clearly imple- 
ment unitary evolutions corresponding to the interaction 
°a ® a p +l - The best way of 'simulating' this interaction 
is hardware-dependent. 

In order to estimate the running time of the algorithm 
one should keep in mind that n(n— 1)/2 interaction terms 



3 Concerning the electrons of solid states, for instance, one is often interested in the energy gaps near the Fermi level W# 
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have to be coupled to the ancilla register (The selection 
of specific coupling does not produce any time overhead 
P,p|JTl|| ) . The time for converting one interaction term 
into a conditional one seems to be independent of n at 
first sight. But here might come in a subtle n-dependence 
as well: it is possible that the simulation of one term 
Hj t k,a,/3 ® o~z has to be more exact if one wants to find 
spectral gaps for large n. Estimating the error in the 
eigenvalues of the total Hamiltonian resulting from an 
error in the simulation of each coupling might be a dif- 
ficult task of perturbation theory jl4j . But the error of 
the total Hamiltonian is only n 2 times the error of the 
single coupling terms. Since we are interested in of gaps 
of the order 1, it should be possible to reproduce them 
with polynomial time overhead. 

The algorithm shows a nice application of 'simula- 
tions of Hamiltonians' Hfj],|||: The algorithm can used to 
find eigenvalues of Hamiltonians by simulating the corre- 
sponding unitary evolution on a quantum computer. 

Apart from its possible application in many-particle 
physics, the algorithm shows a new aspect for under- 
standing the double-role of self-adjoint operators as ob- 
servables and generators of transformations, since we 
have constructed a measurement procedure for H by us- 
ing the evolutions exp(-iHt) as black-box procedures^. 



This gives a new scheme for measuring non-trivial joint 
observables. 

The reader might get confused about one little para- 
dox: we claimed to measure the actual energy value of 
the system despite the fact that only energy differences 
have physical relevance (as long as general relativity is 
not involved). This paradox is resolved by specifying 
more precisely what the algorithm does. We assumed 
the Hamiltonian H to be traceless. In case we start with 
a non-traceless one, only the traceless part will survive 
the conversion procedure H (g) 1 i— ► H <g> a z . Hence the 
algorithm measures the difference between the actual en- 
ergy and the trace of H, i.e., the average of the energy 
values over all the states. This seems again to be para- 
dox: in case the system is in a definite eigenstate, why 
can the ancilla register 'feel' the average energy over all 
the other states despite the fact that the system never 
deviates from its actual eigenstate? This paradox is only 
pretended by the 'infinitesimal language'. Of course the 
system does deviate a little bit from the actual eigenstate 
during the algorithm. Note that the procedure simulat- 
ing H ®o~ z consists of a long sequence of transformations 
(close to the identity) which do not leave the eigenstates 
of H invariant. 
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